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We present the results for Bk calculated using HYP-smeared staggered fermions using lattices 
generated by the MILC collaboration using the asqtad staggered action. We have done the calcula- 
tion on 8 ensembles of these lattices, including three different lattice spacings (a — 0. 12, 0.09, 0.06 
fm). We fit the data to forms based on those predicted by SU(3) mixed-action partially- 
quenched staggered chiral perturbation theory. Our preliminary result is Bk (NDR, pi =2 GeV) = 
0.528 ±0.011 ±0.048, where the first error is statistical and the second systematic. The error 
turns out to be larger than that from an analysis using SU(2) chiral perturbation theory. 
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1. Introduction 

Indirect CP violation violation in the neutral kaon system and is conventionally parameterized 
by e = (2.228 ±0.011) x 10^^*. In the Standard Model, this experimental parameter is related to 
the matrix element of a AS = 2 four-fermion operator, itself parameterized by 

{Ko\sri;dsri;d\Ko) 

(8/3)/|Mi • ^'-'^ 
The relation between experiment and theory takes the form (see, e.g., Ref. [|l|, ^) 



£ = exp(/0e) sin(0e) \/2Ce ImA, X Bk + ^ 
where 0e is the (well-determined) phase of Sk, 

and the dependence on charm and top masses enters through the Inami-Lin functions 50,3 

X = ReAc[T]i5o(xc) - risS3{x,,Xt)] - ReA^T]25o(x,) , x,- = mf/M^ , (1.4) 



+ dim- 8 contributions , (1.2) 



Cs = ^§^, h = V*V,, (1.3) 



with T]i_3 Wilson coefficients. The ^ = ImAo/ReAo term in ( |1.2| ) provides a small correction to the 
proportionality of ex and Bjc, one that must be included now that results for Bfc have errors below 
the 10% level [^. For completeness we note that the contributions of operators of dimension- 8 
have been estimated to be very small, enhancing B/f by ~ 0.5% 

An accurate determination of Bk thus constrains the CKM matrix, as discussed here in the 
review of Van de Water [^. Several groups are pursuing this calculation, using a variety of lattice 
fermions [Q]. Since chiral symmetry constrains the four-fermion matrix element, lattice fermions 
with some measure of chiral symmetry are preferred. Our approach is to use improved staggered 
fermions, which are computationally cheap, although one pays some price due to taste-breaking 
effects. More precisely we use HYP-smeared valence staggered fermions (since HYP-smearing is 
effective at reducing taste-breaking) and asqtad staggered sea quarks (i.e. the MILC ensembles). 
We take advantage of the exact U{1)a chiral symmetry and control the systematic errors due to the 
taste symmetry breaking using staggered chiral perturbation theory (SChPT) ^ |8|]. 

This paper is the first report in a series of the four describing different aspects of our work. 
Here we explain briefly the fitting procedures and results from our SU(3) SChPT analysis. 



2. SU(3) Staggered ChPT Analysis 

In setting up a calculation of Bk using staggered fermions, one must deal with the additional 
taste degree of freedom. For the sea-quarks, this is done by using a rooted determinant. For the 
valence quarks, one must choose a particular taste for the external kaons, and choose the lattice op- 
erator so as to match the matrix element onto the desired continuum one. Away from the continuum 
limit, there are then several sources of error: taste-symmetry breaking coupled with rooting leads 
to unphysical (unitarity violating) effects; taste-symmetry breaking in the valence sector leads to 
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significant changes to the size of chiral logarithms; and approximate matching of the lattice and 
continuum matrix elements (in our case using one-loop perturbation theory) leads to truncation 
errors proportional to powers of a. All these effects are incorporated into SChPT, which thus pro- 
vides a method to systematically control and remove these sources of error. It has been applied 
successfully to the analysis of the n, K, D and B meson masses and decay constants, as reviewed 
inRef. g]. 

The theoretical SChPT analysis was generalized to Bk in Ref. [^], and the result worked out to 
next-to-leading order (NLO). This work assumed the same staggered action for both valence and 
sea quarks and thus is not directly applicable to our mixed-action setup. We have carried out the 



appropriate generalization, as summarized in Ref. []10[]. Full details will be presented in Ref. [|11[]. 

The power-counting used in the SChPT calculation is niq/ Aqcd ~ (a'Aqcd)^ ~ a/TT ~ a^. 
Discretization errors proportional to and to a^a^ are treated as of the same size, since, in practice, 
the latter (which include taste-breaking effects) are numerically enhanced. Since the dominant 
one-loop matching contributions are included, and the coefficients of the taste-breaking operators 
known to be small, these residual one-loop effects (denoted a/n above) are treated as comparable 
to the unknown two-loop effects (a^ above). For more discussion see Ref. 

In the SU(3) SChPT analysis niq can be the mass of the light or strange sea-quarks (mi and 
nts, respectively), or the mass of the down or strange valence quark (nix and niy, respectively). An 
important issue is whether m^/AQCD and my /Aqqu are small enough for NLO ChPT to be useful. 
We address this phenomenologically by including an analytic NNLO term. 

At fixed a, the NLO expression contains 15 different low-energy coefficients (LECs).' 4 of 
these are present in the continuum, while the remaining 1 1 arise from taste- violating discretization 
and truncation errors (including 2 which are introduced by the use of a mixed action). Although 
we calculate Bk for 55 different combinations of m^^ and my, and in some cases for several values 
of mi, it is not possible to do an unconstrained fit to the full functional form. This is partly because 
the functions fall into groups within which they are very similar (differing, for example, by the 
taste of the pion in the loops). The differences within groups are comparable to or smaller than 
our statistical errors. We have chosen to proceed by picking one of functions from each group as 
a representative. In addition, as described in Ref. [|l^], the success of HYP smearing at reducing 
taste breaking impUes that two of the LECs are very small, and so we have dropped these terms. 

After these choices, the number of parameters associated with taste-breaking is reduced from 
11 to 3, so the total number is reduced to the more manageable 7. We write the fitting function as: 



/th = £c,7^- (2.1) 



7 
I 

where c, are the coefficients (related to SChPT LECs), while the functions are 



F2 = G/A\ F, = {G/A^f, F4 = Fj'\2.2) 



i^5 = ^^^, Fe=Fi'l, F,=F^\ (2.3) 



Not including f^, which we fix at its experimental value, taking the scale from the MILC collaboration. 
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and depend on the various pseudo-Goldstone boson (PGB) masses. Our notation is that G = Kp is 
the squared mass of the taste-<§5 valence PGB with flavor composition xy, and 7b are the squared 
masses of the taste-B valence PGBs with compositions xx and yy, respectively, and A), ^ 47ifji is 
the ChPT scale. We set A;^. = 1 GeV in our fits — changing this value results in shifts in the LECs. 

The function Fi contains the continuum-like chiral logarithms, with discretization errors en- 
tering through the taste-breaking in the masses of the PGBs in the loops. The quark-line connected 
part is: 

rB p 

{G + XB)i{XB) + {G + YB)i{YB)+2{G-KB)i{KB)- 2GKJ{Kb ) 



Mconn = L -y 
B=/,Py,AT ^ 



(2.4) 



where is the fractional multiplicity of taste B, and the chiral logarithmic functions are 



m = X\og{X/^l^) , £{X) = = -log(X/M^R) - 1 , (2.5) 

up to (known) finite-volume corrections. The expression for the quark-line disconnected part Mjisc 
is more lengthy: it can be deduced from the results in Ref. [j8|] and will be given in Ref. [pi]]. The 
important point is that, like Mconn it is fully predicted in terms of fji and PGB masses that can be 
calculated in our simulation or obtained from the results of the MILC collaboration. 

The analytic terms F2 and F5 are also present in the continuum, as is the NNLO analytic term 
Ft,. Note that F5 vanishes when nix = niy. 

The remaining three terms are due to taste-breaking discretization or truncation errors, and are 
thus pure lattice artefacts. Of the three, only F4 is non-vanishing when nix = niy. For brevity, we 
show only the two simplest functions: 

Fi'^ = I -^{iGiiKT) + [^(Xt) + ^(Ft) - 2e{KT)] } , (2.6) 
Fit = l^{[i{XT)+i{YT)-2e{KT)]}, (2.7) 
F^^ can be deduced from Ref. |^] and will be given in Ref. [pj|]. 

3. Fitting procedure 

We calculate Bk using standard methods | [T2[ ] on the MILC lattices listed in Table [T[ We use 
the 10 valence quark masses given in Table |^ and thus have 55 different mass combinations for the 
PGBs: 10 degenerate and 45 non-degenerate. We use one-loop matching. 

We use several fitting procedures, but focus here on the one used for our central value. In 
this fit we do not yet include the finite-volume corrections predicted by ChPT. First, we fit the 



degenerate data. The general form (2J_) reduces to ^^^j c,F,- in the degenerate limit. We know the 



expected size of each term; in particular, we expect C4 to have the magnitude C4 a^Aq^^) if it is 



dominated by discretization errors. We implement this prior knowledge using Bayesian fitting [ |13| ] 
in which the standard is augmented by the addition of 

Xprior ^2 ' (3-1) 

t>4 
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a (fm) 


ami /arris 


geometry 


ens X meas 


ID 


0.12 


0.03/0.05 


20^ X 64 


564 X 1 


CI 


0.12 


0.02/0.05 


20^ X 64 


486 X 1 


C2 


0.12 


0.01/0.05 


20^ X 64 


671 X 9 


C3 


0.12 


0.01/0.05 


28^ X 64 


274x8 


C3-2 


0.12 


0.007/0.05 


20^ X 64 


651 X 1 


C4 


0.12 


0.005/0.05 


24^ X 64 


509 X 1 


C5 


0.09 


0.0062/0.031 


28^ X 96 


995 X 1 


Fl 


0.06 


0.0036/0.018 


48^ X 144 


513 X 1 


SI 



Table 1: MILC lattices used for the numerical study. Here "ens" and "meas" are the number of gauge 
configurations and measurements per configuration, respectively, and "ID" is a label. 



a (fm) 


arrix and arriy 


n 




0.12 


0.005 X n 


1 2 3 


.,10 


0.09 


0.003 X n 


1,2,3,.. 


.,10 


0.06 


0.0018 X n 


1,2,3,.. 


.,10 



Table 2: Valence quark masses used for the numerical study (with nix < niy). 

1 — I — I — r 



N-BT7 



S (D ffl 

■(^■■■^■■(jj-i-§-g-g-g-g-@-@-i-g-@ 

o 



I I I I I I I I I I I I I I I I I I I I I I I I 
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m^'^ (deV^) (OeV^) 

Figure 1: Bk vs. m|. (left) and A.Bk vs. m\ (right): The results are on the C3 lattices. The fit types are 
D-BT7 (left, degenerate case) and N-BT7 (right, non-degenerate case). 




We take a^ =0 since we do not know the sign of C4, and 04 = a'^K^^^ with Aqcd = 0.3GeV. 
Minimizing the augmented leads to the fit we label D-BT7 (with "D" for degenerate).^ An 
example is shown in the left panel of Fig. |]. 

Next we fit /th to all 55 data points, but using the results for ci_4 and their errors from the 
D-BT7 fit as priors for the first four parameters. We also include priors for the coefficients C(, and 

^At the present stage, our includes only the diagonal elements of the correlation matrix. 
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cj, whose expected magnitudes are cg ~ C4 and c^ w Aq(2dC4- In total then, is augmented by 

^''~f (3.2) 

i=l i=6,7 

with ai_4 set to the central values of ci_4 from the D-BT7 fit, and d'i_4 set to the corresponding 
errors, while = a? = 0, dg = o^Aqq^ and dj =a ^A^^p. The resulting fit is labeled N-BT7 ("N" 
for non-degenerate). An example is shown in the right panel of Fig. |l|, where ABk is the residual: 

ABj,=Bjf(data)-/,h. (3.3) 

The fit works well except for the non-degenerate points in which rrix takes its smallest value, which 
we suspect is due in part due to finite-volume effects. 

We have tried many variations of this fitting strategy: relaxing the constraints on 04.6,7 to the 
values one might expect if O(a^) [rather than 0{a^)] effects were dominant (type N-BT7-2); fitting 
in one stage to all data-points with either 0{a^) or 0{a^) priors on 04,5,7 (types N-BT8 and N-BT8- 
2); and simply fitting without priors (N-T2). The results for the continuum-like parameters ci_3 
and C5 are consistent between fits, while those for the lattice artefacts C4, C6-7 vary significantly. 
An intriguing result is that we find cj to be remarkably smaller than the expected size [C5 0(1)]. 

For each fit, we determine a value for Bk at the physical valence-quark masses by setting 
all lattice artefacts to zero in the SChPT expression. We can then attempt to extrapolate to the 
continuum and infinite volume limits, and to the physical values of mt and m^. Examples of some 
of these extrapolations are shown in the companion proceedings [ 14, p3|]. 



4. Error budget and conclusion 



source 


error (%) 


description 


statistics 


2.0 


N-BT7 fit 


discretization 


2.5 


diff. of (SI) and (a = 0) 


fitting (1) 


0.6 


diff. of N-BT7 and N-BT8 (C3) 


fitting (2) 


5.4 


diff. of N-BT7 and N-BT7-2 (C3) 


fitting (3) 


1.3 


fit w.r.t. ami 


finite volume 


2.2 


diff. of C3 and C3-2 


matching factor 


6.3 


Ab'^k^' (SI) 


scale ri 


0.07 


uncertainty in ri 



Table 3: Preliminary error budget for Bk obtained using SU(3) SChPT fitting. 



In Table ^ we summarize our present best estimates of the uncertainties in Bk coming from 
various sources. These are conservative estimates, some of which may decrease with further anal- 
ysis of our present data set. The method by which we estimate these errors is outlined in the 
"description". The first fitting error estimates the (rather small) impact of fitting in two stages 
rather than in one, while the second fitting error estimates the impact of changing the widths 04,6,7 
for the priors. The latter error is one of our two dominant systematic errors — it indicates that we 
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cannot pin down the functional form of Bk well enough to extrapolate to the physical light-quark 
mass with less than 5% uncertainty (at least on the coarse lattices). The other dominant error is that 
due to the truncation of the matching factor at one-loop order. This error is discussed in one of the 
companion proceedings [p^]. 

Our current, preliminary estimate of Bk using SU(3) ChPT fitting is 

B^(NDR,Ai = 2 GeV) = 0.528 ±0.011 ±0.048 [SU(3), PRELIMINARY.] (4.1) 

where the first error is statistical and the second combines all systematic errors. This turns out to 



be less accurate than the results from SU(2) ChPT fitting [14]. 
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